Noname manuscript No. 

(will be inserted by the editor) 



P. Giommi • S. Colafrancesco 



Non-thermal Cosmic Backgrounds 
and prospects for future high-energy 
q ■ observations of blazars 



O 

<N 

Received: date / Accepted: date 



Abstract We discuss the contribution of the blazar population to the ex- 
tragalactic background radiation across the electromagnetic (e.m.) spectrum 
with particular reference to the microwave, hard-X-ray and 7-ray bands. 
Our estimates are based on a recently derived blazar radio LogN-LogS that 
was built by combining several radio and multi-frequency surveys. We show 
that blazar emission integrated over cosmic time gives rise to a considerable 
broad-band non-thermal cosmic background that dominates the extragalac- 
tic brightness in the high-energy part of the e.m. spectrum. We also estimate 
the number of blazars that are expected to be detected by future planned or 
hypothetical missions operating in the X-ray and 7-ray energy bands. 

■ 

Keywords galaxies: active - galaxies: ■ blazar: BL Lacertae surveys 

6 
& 



(N 
(N 
O 



1 Introduction 



Active Galactic Nuclei (AGN) are well known to dominate the high-energy 
(soft-X-ray and beyond) high Galactic latitude sky. Their radiation inte- 
grated over cosmic time can explain most, if not all, the extragalactic back- 
ground radiation (e.g. [4Tl24U17ll7| ). Historically, AGN have been classified in 
many, sometime inconsistent or confusing, ways, depending on how they ap- 
peared in surveys performed in different energy bands and flux limits, often 
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determined by the technology available at the time of the survey, or depend- 
ing on some observational parameter like e.g. the flux ratio between radio 
and optical emission, the equivalent width of their emission lines, etc. 

In the following we use the widely accepted standard paradigm where 
AGN are powered by accretion onto a super-massive black hole, e.g. [301311 
EH], to divide this class of sources, into two broad categories defined by their 
intrinsic emission mechanism. 

AGN whose power is dominated by radiation of thermal origin 

produced by matter that is strongly heated in the inner parts of a disk 
of matter falling onto the super-massive black hole and illuminates the 
circumnuclear matter that is responsible for the broad line emission. We 
call these objets Thermal Emission Dominated AGN or TED-AGN. 
AGN whose power is dominated by Non-Thermal radiation (or 
NTED- AGN) where most of the emission is generated through non-thermal 
processes like the synchrotron and the inverse Compton mechanism by 
particles accelerated in a jet of material that moves at relativistic speed 
away from the central black hole. The jet itself if formed converting part 
of the accretion energy in a way that is presently not well understood. 

Within this general definition the class of blazars corresponds to the small 
subset of NTED- AGN that are viewed at a small angle w.r.t. the jet axis (for 
this reason their emission is strongly amplified by relativistic effects |2"ll30| L 
whereas Radio Galaxies are those NTED-AGN that are viewed at a large 
angle w.r.t. the jet axis. Here we do not distinguish between line-less objects, 
or BL Lacs, and broad-line Flat-Spectrum Radio Quasars, or FSRQs. 

When the accretion and jet emission coexist in the same object and pro- 
duce similar amounts of radiation we have AGN that are a mixture of the 
TED and NTED type that appear like broad- line radio galaxies (e.g. 3cl20). 

The overall cosmic background energy has two well understood compo- 
nents: the primordial black body emission peaking at microwave frequencies, 
or CMB, and the X-ray apparently diffuse emission arising from the accre- 
tion onto super-massive black holes in AGN integrated over cosmic time, or 
CXB. We will show that blazars add a third non-thermal component that 
at low frequencies contaminates the CMB fluctuation spectrum and compli- 
cates its interpretation ^5 , while at the opposite end of the e.m. spectrum it 
dominates the extragalactic background radiation |7j. 

In the following we estimate the contribution of blazars to non-thermal 
cosmic backgrounds starting from a recently derived deep radio LogN-LogS. 
The broad-band electromagnetic spectrum of a blazar is composed of a syn- 
chrotron low-energy component that peaks [in a Log(v) — Log(v f (v)) rep- 
resentation] between the far infrared and the X-ray band, followed by an 
Inverse Compton component that has its maximum in the hard X-ray band 
or at higher energies, depending on the location of the synchrotron peak, and 
extends into the 7-ray or even the TeV band. Those blazars where the syn- 
chrotron peak is located at low energy are usually called Low energy peaked 
blazars or LBL, while those where the synchrotron component reaches the 
X-ray band are called High energy peaked blazars or HBL . LBL sources 
are the large majority among blazars |2()| and are usually discovered in ra- 
dio surveys, while HBL objects are preferentially found in X-ray flux limited 
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surveys, since at these frequencies they are hundreds, or even thousands, of 
times brighter than LBLs of similar radio power. 



2 The radio LogN-LogS of blazars and their contribution to 
extragalactic cosmic backgrounds 

The deep blazar radio (5GHz) LogN-LogS shown in Fig. ^ has been recently 
derived by [J] who showed that it can be described by a broken power law 
with parameters defined in eq. l|T]): 



□ 2Jy Flat Spectrum Radio Survey 
O Uy ASDC-RASS-NVSS 
★ WMAP selected Blazars 
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Fig. 1 The blazars radio (5GHz) LogN-LogS of 7 who built it combining several 
surveys as indicated in the top right part of the plot (see |7j for details). 



M( . Q ,_ J 5.95 1CT 3 x S- 1 - 62 S> 0.015 Jy m 
iV ^>^-\ O.f 25 x S~ - 9 S < 0.015 Jy W 

Once the LogN-LogS of a population of sources is known in a given energy 
band it is possible to estimate their emission in other parts of the electro- 
magnetic spectrum, provided that the flux ratio between the two bands is 
known. In this section we deal with flux ratios and Spectral Energy Distri- 
butions (SED) of blazars that will be used to estimate the contribution to 
cosmic backgrounds at microwave frequencies and in the X-ray and 7-ray 
bands. 

The distribution of f x /f T flux ratios shown in Fig|2] (see [7] for details) 
spans about four orders of magnitudes implying that the X-ray flux of blazars 
with the same radio flux can be different by up to a factor 10,000! 

In the following we estimate the blazar contribution to the cosmic back- 
grounds basing our calculation on the radio LogN-LogS of Fig. ^ and on flux 
ratios in different bands (Fig [21 or on observed blazar SED. 
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Fig. 2 The / x //r distribution of blazars estimated by 
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Fig. 3 The extragalactic Cosmic Background Energy distribution at microwave, 
X-ray and 7-ray energies. Th e CMB is represented by a black-body spectrum with 
temperature of 2.725 °K 1161 ; the CXB is taken from HEAO-1 measurements 1151 
|3| and has been scaled to match the more recen t BeppoSAX, ASCA, and XMM- 
Newton results in the 2-10 keV band I32I14I12I . The gam ma ray background is 
derived from the COMPTEL data in the range 0.8-30 MeV [TT] and from EGRET 
data in the range 30 MeV - 50 GeV. We rep ort the results of two d iffer ent analyses 
of the EGRET data: open circles from |26| and filled circles from |27|. which uses 
an improved model of the Galactic diffuse continuum gamma-rays. As for the TeV 
diffuse background, we report the upper limits in the 20 - 100 TeV region derived 
from the HEGRA air shower data analysis 0. In the 1 TeV - 1 PeV energy range, 
other experiments give only upper limits and there is no clear observation of a 
diffuse photon signal yet (see [7j for more detail). 



2.1 The Cosmic Microwave Background 



The contribution of blazars to the CMB has been estimated in the past from 
different viewpoints (e.g. 29,5 ). Here we calculate the integrated microwave 
background intensity as follows 

r 1J v dN 

hlazars = S — — dS , (2) 

JO.lmJy "<J 
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where dN/dS is the differential of eq. The minimum integration flux of 0.1 
mjy for S m i n is likely to be conservative, since blazars with radio flux near or 
below 1 mjy are already included in the Einstein Medium Sensitivity Survey 
BL Lac sample |23| . The integrated intensity Iuazars is converted from 5GHz 
to microwave frequencies by convolving the flux value with the observed 
distribution of spectral slopes between 5GHz to microwave frequencies, as 
estimated from the Uy-ARN sample (see jS] for details). 




Fig. 4 The contribution of blazars to the CMB fluctuation spectrum in the WMAP 
41GHz (left panel) and 94 GHz (right panel) channels, as evaluated from the LogN- 
LogS given in Fig. (solid line) . We also show the angular power spectrum for 
the blazar population by adding an estimate of the possible contribution of radio 
sources with steep-spectrum at low radio frequencies which flatten at higher fre- 
quencies (dashed line) . The dotted line also includes the effect of spectral and flux 
variability. Although this additional contamination may be substantial a precise 
estimation can only be done through simultaneous high resolution observations at 
the same frequency. A typical CMB power spectrum evaluated in a TCDM cos- 
mology with Q m — 0.3, Qa = 0.65, Qb = 0.05 which best fits the available data is 
shown for comparison. 



The fractional contamination, Iuazars/ Icmb of blazars to the CMB in- 
tensity {Icmb) and the corresponding apparent temperature increase at fre- 
quency of 90-300 GHz has been derived by [7] . Blazars also contribute to the 
temperature anisotropy spectrum of the CMB with a spectrum 

(AT) blazar = [^^Cdil+l)] 1 ' 2 , where 

C^blazar = J dS —j^ S 2 , (3) 

The quantity on the right hand side in eq. J3J) is the usual Poisson shot-noise 
term, |21U28j and we neglect here the clustering term, w(/) 2 , which adds to 
the Poissonian term since there is not a clear estimate of the blazar clustering 
yet. Note, however, that the inclusion of clustering significantly increases the 
amount of contamination, especially at large angular scales |S]. 

For the blazar population described by the LogN-LogS given in Fig. ^ 
we found CV, blazar ~ 54.4 Jy 2 sr _1 at 41 GHz and C^biazar ~ 51.6 Jy 2 sr _1 at 
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94 GHz. When translated into temperature units (see [2]) we found a value 
Cf,biazar « 2.22 • 1CT 2 /iK 2 sr at 41 GHz and C^biazar « 1.09 • 10~ 3 ^K 2 sr at 
94 GHz. We show in Fig. 0] the quantity (zAT)biazar and compare it to the 
CMB fluctuation power spectrum which best fits the available data. These 
values (shown as solid lines) have to be considered as a definite lower limit 
for C^biazar, since they neglect the contribution of steep-spectrum sources at 
low frequencies which flatten at these frequencies (41 and 94 GHz, dashed 
lines) and the effect of flux variability (dotted lines). 

The blazar flux variability at millimeter wavelengths may be substantial 
(higher than factors 3-10 on time scales of weeks to years seen at cm wave- 
lengths) and could increase the amount of contamination of CMB maps when 
these are built over long integration periods. 

The contamination level shown in Figs. 01 is the one expected in the case 
where no blazars are removed from the CMB data. We expect that a correct 
procedure to derive the CMB power spectrum that takes into account the 
full point-like source contribution implied by our LogN-LogS, would both in- 
fluence the shape of the expected power spectrum and increase the statistical 
uncertainties of the WMAP data points, especially at high multipoles, where 
the blazar contribution is larger. 

The previous calculations we performed neglecting the blazar clustering 
and thus they must again be considered as a lower limit to the realistic angu- 
lar power spectrum contributed by these sources. The effects of clustering on 
the CMB fluctuation spectrum has been partially estimated by some authors: 
(e.g. EEi}). Expectations for the clustering effect strongly depend both on 
the adopted model for the source counts and on their clustering model. Based 
on the correlation function of ^Sj on that for the SCUBA sources and 
on our blazar LogN-LogS, we expect that the contamination of the first peak 
of the fluctuation spectrum (at the WMAP 41 GHz channel) is at a level in 
the range 20-25 %. This estimate does not include possible variability effects 
and additional core-dominated radio sources 



2.2 The Soft X-ray background 

We estimate the contribution of blazars to the CXB at 1 keV using two 
methods: i) converting the integrated radio flux calculated with Eq. J2J into 
X-ray flux with the observed distribution of / x // r flux ratios in Fig. |2 and 
ii) converting the integrated blazar contribution to the CMB (at 94 GHz) to 
X-ray flux using the distribution of a^ x , the microwave (94 GHz) to X-ray 
(IkeV) spectral slope, estimated from the subsample of blazars detected by 
WMAP in the 94 GHz channel, for which an X-ray measurement is available 
(see Fig. 9 in [7]). The first method gives a total blazar contribution to the X- 
ray background of 2.7 x 10~ 12 erg cm _2 s _1 deg _1 (about 70% of which is due 
to HBL sources with / x // r > 5 x 10" 12 erg cm^s" 1 Jy" 1 ) in the ROSAT 0.1- 
2.4 keV energy band. Assuming |7j an average blazar X-ray energy spectral 
index of a x = 0.7, this flux converts to 2.6 x 10~ 12 erg cm _2 s _1 deg _1 in 
the 2-10 keV band or 11% of the CXB, which is estimated to be 2.3 x 10 -11 
erg cm _2 s _1 deg _1 \'2'2\ . 
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The distribution of a^x has an average ( a^ x ) — 1.07 and a standard 
deviation of 0.08 corresponding only to about a factor 3 in flux. This dis- 
tribution is much narrower than the one between the radio and X-ray band 
(Fig. while the dispersion is comparable to that expected from blazar 
variability at radio and especially at X-ray frequencies. The X-ray flux can 
therefore be estimated simply as /ifc e V=l-4x 10~ 7 f^GHz (see, e.g., jjj). 

Since the blazar integrated emission at 94 GHz is 7.2 x 10~ 6 CMBg4G£f z or 
0.64 Jy/deg 2 and the cosmic X-ray background is 2.3 x 10 -11 erg cm _2 s _1 [22] 
in the 2-10 keV band (equivalent to 2.31/iJy/deg 2 at 1 keV), we have that 
i lkeV = 0.09 ^Jy/deg 2 or 3.9% of the CXB for / 94 G^=0.64Jy/deg 2 . Con- 
sidering that the a^x distribution of 7 only includes LBL objects and that 
HBL sources make up two thirds the total contribution to the CXB, the per- 
centage scales to about 12% which is very close to the 11% obtained with the 
previous method. Both results are in good agreement with the independent 
estimate of who find that the radio loud AGN content of the CXB is 13% 
in the XMM-Newton bright serendipitous survey. 

2.3 Hard X-ray - soft 7-ray Background 

The number of sources detected at energies larger than soft X-rays is still 
rather low, so building reliable distributions of flux ratios between radio or 
microwaves and the Hard X-ray/7-ray fluxes is not currently possible. Thus, 
in order to estimate the blazar contribution to high energy Cosmic Back- 
grounds (E > 100 keV), we therefore followed a different approach: we 
extrapolated the predicted blazar integrated intensity at microwave frequen- 
cies (eq. El to the hard X-ray and soft 7-ray band using a set of hypothetical 
SSC spectral energy distributions. 

Figure0(left panel) shows the CMB, CXB, and CGB together with three 
predicted SEDs from a simple homogeneous SSC models whose parameters 
are constrained to 1) be consistent with the expected integrated flux at 94 
GHz, 2) have the a^ x slope equal to the mean value of the WMAP blazars 
{a^x = 1-07), and 3) possess a radio spectral slope equal to the average value 
in the WMAP sample. The three curves, so constrained, are characterized by 
synchrotron peak frequencies of v pea k = 10 12 8 , 10 13 5 , and 10 138 Hz. A high 
value of v pe ak largely overestimates the observed hard- X-ray to soft 7-ray 
(« 10 20 — 2 x 10 22 Hz or w 500 keV-10 MeV) cosmic background, whereas 
a too low value of u pea k predicts a negligible contribution (see Fig. left 
panel). The case v pea k — 10 13 5 Hz predicts 100% of the cosmic background. 
Since the Log(^ pea fc ) values of blazars in the WMAP and other catalogs peak 
near 13.5 and range from 12.8-13.7 within one sigma from the mean value, 
we conclude that blazars may be responsible for a large fraction (possibly 
100%) of the hard-Xray/soft 7-ray cosmic background. 

2.4 7-ray Background 

The SSC distributions of Fig. predict a negligible blazar contribution to 
the extragalactic 7-ray Background above 100 MeV. Nevertheless, it is well 
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known that blazars are the large majority of the extragalactic 7-ray (E > 100 
MeV) identified sources detected by the EGRET experiment therefore 
they are likely to contribute to the 7-ray background in a significant way. 
Indeed, JS| concluded that blazars should make a large fraction, if not the 
totality, of the extragalactic 7-ray background. However, these early calcula- 
tions relied upon a very small sample and had to assume no strong variability, 
a characteristic that was later demonstrated to be extremely common in 7- 
ray detected blazars. 

Figure 11 in 7 compares the energy distribution of the CMB, CXB, 
and CGB to the SED of 3C 279, a well-known bright blazar detected by 
EGRET, scaled to the integrated blazar flux intensity as calculated with 
eq. © ■ Taking into account the strong variability of blazars seen across the 
whole e.m. spectrum jj] show that while the contribution to the CXB can 
range from a few percent to over 10% in the higher states, the predicted flux 
at 7-ray frequencies ranges from about 100% to several times the observed 
cosmic background intensity. Such large excess implies that either sources 
like 3C279 are not representative of the class of blazars or their duty cycles 
at 7-ray frequencies are very low. A way of quantifying the ratio between 
the 7-ray intensity predicted by assuming that the source is representative 
of the entire population and the actual background intensity is to use the 
microwave (94 GHz) to 7-ray (100 MeV) slope a M7 = -gjte^M^. 




Log frequency > (Hz) Log freqtencj if (Hz) 



Fig. 5 Left: The possible contribution of LBL blazars to the cosmic background. 
Simple SSC energy distributions fail to reproduce the observed intensity and slope 
of the 7-ray background. The observed range of 7-ray emission with EGRET, nor- 
malized to radio flux, indicates that the duty cycle at these energies must be low 
(see text for details). Right: The spectral slopes between the 5GHz emission and 
the 100 keV and 1 MeV high energy fluxes in the SSC distribution of LBL blazars 
that reproduces the observed X and 7-ray cosmic background. 

For comparison, we defined ^ 7ioo%CGB = 0.994 as the a M7 of a hypo- 
thetical source that would produce 100% of the CGB if representative of the 
entire class. Any real source with a M7 flatter than 01^ — 0.994 would 

then predict an integrated flux in excess of the observed 7-ray background if 
representative of the entire population; alternatively, its duty cycle must be 
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lower than 100%. Fig. (left) shows the minimum and maximum values of 
a M7 among the sources detected by the EGRET experiment. 

2.5 TeV background 

All blazars detected so far at TeV energies are of the HBL type, therefore 
we estimate the blazar contribution to the TeV background as in the case 
of the soft 7-ray Background but only considering the HBL component and 
not the entire blazar population. These extreme objects have been estimated 
to be about 0.1 % of the blazar population [7|. Using the SED of the well 
known TeV blazar MKN421 normalized at 94 GHz (see Fig.13 in 0), so that 
the flux is scaled to 0.1% of the intensity produced by the entire population 
of blazars, we see that, despite HBLs are a tiny minority, their integrated 
X-ray flux makes up a fairly large fraction of the CXB and that their TeV 
emission may produce a significant amount of extragalactic light. We note, 
however, that since extreme HBLs, such as those of the Sedentary survey, are 
very rare (one object in several thousand degrees with flux above a few mJy), 
the extragalactic light at TeV energies should be very patchy, associated to 
single sources, rather than a diffuse light resulting from the superposition of 
many unresolved discrete sources. 

3 Prospects for future high-energy observations of blazars 

We have seen that blazars, and in general NTED-AGN, emit non-thermal 
radiation across the entire electromagnetic spectrum. A clear understanding 
of the blazar content of the extragalactic high-energy sky is necessary as new 
missions are planned or are studied by European and other space agencies to 
explore in more detail the hard X-ray and 7-ray bands and to open the last 
energy windows (e.g. the few hundred keV-few MeV region) that for technical 
difficulties have so far remained almost unexplored but that are extremely 
important for the understanding of the physical processes powering NTED- 
AGN. 

In the following we apply the same method used in the previous sections to 
estimate the number of blazars that could be detected by future microwave, 
X-ray and 7-ray observatories. We report the results in Tablc|3whcrc column 
1 gives the energy band considered, column 2 gives the limiting sensitivity 
reached in that energy band, column 3 gives the corresponding 5GHz flux, 
and column 4 gives the number of blazars expected in the high Galactic 
latitude sky (« 25,000de.g 2 ). 

Figure [2] shows that the distribution of / x / / r ratios is extremely broad 
reflecting the fact that the flux in the X-ray region is deeply affected by 
the position of the synchrotron peak in Log{y) — Log(y j (f)) plots. Blazars 
where the peak is located in the far infrared (LBL) have very little X-ray flux 
whereas extreme HBL objects can be thousands of times brighter at X-ray 
frequencies than LBLs of the same radio flux. Figure [21 implies that a survey 
as deep as 10 -15 erg cm _2 s _1 would detect the whole / x //r distribution 
of blazars with 5GHz flux larger than 10 mJy, corresponding to well over 
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Table 1 The number of blazars in the high Galactic latitude sky at different 
sensitivities and energy bands 



Energy band 


Flux limit 


Flux limit 


No. of sources 


(energy band) 


(5GHz) 


at high b| 








(25,000 deg 2 ) 


Micro- wave 


70 mJy 


100 mJy 


5000 


(100 GHz) 


7 mJy 


10 mJy 


>100,000 


Soft X-Ray 


10" ib erg cm" 2 s _1 


~0.1-10 mJy 


>100,000 


(0.1-2.4 keV) 








Hard X-Ray 


5 W 1 ' 2 erg cm'^s" 1 


1000 mJy 


130 


50 - 150 keV 


1 10~ 13 erg cm" 2 s _1 


20 mJy 


^70,000 


Soft 7-ray 


5 10~ v ph cm-'s" 1 


1000 mJy 


130 


3 - lOMeV 


1 10 -7 ph cm^s" 1 


200 mJy 


2000 


7-ray 


1 10" Y ph cm- z s _1 


Monte Carlo 


w 100 


100 MeV-100 GeV 


3 10" 9 ph cm-V 1 


simulation 


» 3000 



100,000 sources in the high Galactic latitude sky. From Fig. [5] (right panel) 
we can see that the average slope between the 5GHz radio flux and the hard- 
X-ray band and the soft 7-ray band is 0.85 and 0.84 respectively. Since the 
ratio between the radio flux and these energy bands is not expected to have 
a large spread around the mean value like / x // r , the average slopes can be 
directly converted to predict the hard- X-ray and the soft 7-ray fluxes. We 
have reported these predictions in column 2 of Table |31 
A deep all sky survey in the soft X-ray band: 0.1-10 keV. From 
Table we see that a deep all-sky X-ray survey in the soft X-ray band would 
allow us to build a sample of > 100,000 blazars that are expected to be 
above the limiting sensitivity of PLANCK. Given the significant impact of 
the blazar foreground emission on the CMB power spectrum it is important 
to remove this contaminating component from the CMB as much as possi- 
ble. An efficient way to achieve this is to exploit the fact that the spectral 
slope distribution between microwave and soft X-ray flux of LBLs is very 
narrow (see Fig. 9 of [7]) with a dispersion that is probably mostly due to 
intrinsic variability. The soft X-ray flux of LBLs (that is >90% of the blazar 
population) is therefore a very good estimator of the flux at microwave fre- 
quencies and could be used to locate and remove foreground blazars from the 
PLANCK and other CMB missions. 

The extremely large sample of blazars produced by such a survey could 
also be used to study the statistical properties of blazars in great detail, 
including the spatial correlation function, and would identify a very large 
number of targets for the next generation of 7-ray observatories such as AG- 
ILE, GLAST, and future instruments operating in the still poorly explored 
MeV spectral region. 

The expected number of blazars at 7-ray energies above 100 MeV cannot 
be estimated with a simple extrapolation of the 5GHz flux (see Sect. 2.4) and 
therefore we have calculated the preliminary expectations listed in the last 
row of Table |31 through a Monte Carlo method that makes use of the radio lu- 
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minosity distribution of blazars and predicts 7-ray fluxes taking into account 
the a A17 distribution of EGRET detected blazars. This method is currently 
under development as part of our contribution to the GLAST project. 

Acknowledgements It is a pleasure to thank the Organizing and Advisory Com- 
mittee of this meeting for the invitation to present our results in the context of the 
planned and future gamma-ray experiments. 
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